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Abstract 

Modified gravity and generalized second law (GSL) of thermodynamics are interesting 
topics in the modern cosmology. In this regard, we investigate the GSL of gravitational 
thermodynamics in the framework of modified Gauss-Bonnet gravity or /(G)-gravity. We 
consider a spatially FRW universe filled with the matter and radiation enclosed by the dy¬ 
namical apparent horizon with the Hawking temperature. For two viable /(G) models, we 
first numerically solve the set of differential equations governing the dynamics of /(G)-gravity. 
Then, we obtain the evolutions of the Hubble parameter, the Gauss-Bonnet curvature invari¬ 
ant term, the density and equation of state parameters as well as the deceleration parameter. 
In addition, we check the energy conditions for both models and finally examine the validity 
of the GSL. For the selected /(G) models, we conclude that both models have a stable de Sit¬ 
ter attractor. The equation of state parameters behave quite similar to those of the ACDM 
model in the radiation/matter dominated epochs, then they enter the phantom region be¬ 
fore reaching the de Sitter attractor with u> = — 1. The deceleration parameter starts from 
the radiation/matter dominated eras, then transits from a cosmic deceleration to accelera¬ 
tion and finally approaches a de Sitter regime at late times, as expected. Furthermore, the 
GSL is respected for both models during the standard radiation/matter dominated epochs. 
Thereafter when the universe becomes accelerating, the GSL is violated in some ranges of 
scale factor. At late times, the evolution of the GSL predicts an adiabatic behavior for the 
accelerated expansion of the universe. 
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1 Introduction 


Various cosmological observations, coming from the type la supernovae (SNela) surveys [1], the 
large scale structure (LSS) [2], the cosmic microwave background (CMB) anisotropy spectrum 
[3, 4] and the Hubble parameter H(z ) [5], have indicated that the universe is in a phase of 
accelerated expansion. Regarding the accelerated expansion of the universe, there are two main 
categories of probable solutions. One is to assume that in the context of general relativity 
(GR), the universe is dominated by a new cosmic fluid with negative pressure. This kind of 
exotic matter which violates the strong energy condition is so called “dark energy “ (DE) [6, 7]. 
Another alternative, originates from the modification of gravity [8]. In modified gravity (MG) 
theories, there is no require for DE with exclusive properties, but instead, the action contains 
a general function of invariants obtained from the Riemann curvature tensor such as the Ricci 
scalar, R, [9] or the Gauss-Bonnet invariant term, G, [10] or the torsion scalar, T [11], Moreover, 
in [12] it was shown that MG may serve as dark matter (DM). 

One of interesting alternative theories of gravity is modified Gauss-Bonnet gravity, so-called 
f(G)- gravity, where f(G ) is a general function of the Gauss-Bonnet curvature invariant term 
G = R 2 — 4 R^R^ + R.f_ l up f jR ,lvp,J [10, 13, 14, 15, 16, 17, 18]. The f(G )~gravity can justify 
the present accelerated expansion of the universe without resorting to DE. Besides, it can also 
describe the phantom divide line crossing as well as the cosmic transition from deceleration 
to acceleration phase [19, 20, 21]. The cosmologically viable f(G) models need to be close to 
the ACDM model in the deep matter era, but the deviation from it becomes important at late 
times on cosmological scales. An appreciable deviation from the ACDM cosmology yields the 
modification of the matter power spectrum, which can be used as a crucial tool to distinguish 
f(G )-gravity from the ACDM model. 

Additionally, the thermodynamical interpretation of gravity is one of another interesting top¬ 
ics in modern cosmology. From the viewpoint of the physics of the black holes, there is a deep 
rooted connection between thermodynamics and gravity [22]. This connection was first discov¬ 
ered in the Einstein gravity for the Rindler spacetime [23]. It was also shown that by assuming 
the geometric entropy given by a quarter of the apparent horizon area of a Friedmann-Robertson- 
Walker (FRW) universe, the Friedmann equation in the Einstein gravity can be written in the 
form of the first law of thermodynamics [24], The connection between thermodynamics and 
gravity has also been investigated in /(R)-gravity and scalar-tensor theory [25], f(T )-gravity 
[26], Lovelock theory [27] and braneworld scenarios (such as DGP, RSI and RSII) [28]. 

It is also of great interest to generalize our discussion to study the generalized second law 
(GSL) of thermodynamics. The GSL states that the entropy of matter inside the horizon of the 
universe plus the geometric entropy of the horizon is non-decreasing with time [24], Note that 
the ordinary second law of thermodynamics only deals with the entropy of matter inside the 
universe. In the Einstein gravity, it was shown that the GSL in the presence of DE is always 
satisfied [29]. The validity of the GSL was also examined in different theories of MG [29]-[37]. 
Here our main aim is to explore the GSL and the thermodynamics of the apparent horizon in 
f(G)~ gravity, and obtain the condition for the GSL to be satisfied. The paper is structured as 
follows. In section 2, we briefly review the f(G)- gravity. In section 3, we investigate the GSL of 
thermodynamics on the dynamical apparent horizon with the Hawking temperature. In section 
4, we study the dynamics of f(G)- gravity. In section 5, we give numerical results obtained for 
the evolution of some cosmological parameters, check the energy conditions and examine the 
validity of the GSL for two viable f(G) models. Section 6 is devoted to conclusions. 


2 



2 THE f(G ) THEORY OF GRAVITY 

The action of modified Gauss-Bonnet gravity is given by [14]: 

I = J d ^XyJ—g ^7^2 R + / (^) + L r + , ( 1 ) 

where k 2 = 8nG n = 1 and G n is the Newtonian gravitational constant. Also g , R , L r , L m 
and /(G) are the determinant of metric , Ricci scalar, the matter Lagrangian, the radiation 
Lagrangian and a general function of the Gauss-Bonnet term, respectively. The Gauss-Bonnet 
curvature invariant term is defined as 


G = R 2 — 4 + R^R^. 


(2) 


Taking variation of the action (1) with respect to g^ leads to the field equations 


Tpy — Rpu 2 ^^^ [77/ipt/<r T R/nyfJrrii Rpagufi Rpu9ap T Rpagvp 
+ ^R{9pu9ap - 9pa9vp)]y p y a /g + (G/g - f)9pv, 


(3) 


where /g = d//dG. Also Tp V and V p are the energy-momentum tensor and the covariant 
derivative operator, respectively. Now we consider a spatially flat universe described by the 
FRW metric 


ds 2 = -dt 2 + a 2 (t) (dx 2 + d y 2 + dz 2 ') , (4) 

where a(t) is the scale factor. Consequently, we have 

R = 6 (H + 2FT 2 ) , G = 24 H 2 (h + H 2 ) , (5) 

where H = a/a is the Hubble parameter and an overdot stands for a derivative with respect 
to the cosmic time t. In terms of the deceleration parameter q = —1 — H/H 2 , Eq. (5) can be 
written as 

R = QH 2 (1 - q) , G = 24 H*q. (6) 

Substituting the FRW metric (4) into the field equations (3) yields the Friedmann equations 
in /(G)-gravity as 

3 H 2 = G/g - / - 24H 3 /g + Pr + Pm, (7) 


-2 H = -8H 3 / g + 16HHf G + 8 H 2 f G + | p r + p m , (8) 

where p m and p v are the energy density of matter and radiation, respectively. 

The Friedmann equations (7) and (8) can be rewritten in the standard form as [38] 

H 2 = \pu (9) 

H = ~{pt+p t ), (10) 
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where p t and p t are the total energy density and pressure defined as 


Pt = Pm + Pi- + Pg, (11) 

Pt=Pm+Pr+PG • (12) 

Here p G and p G are the energy density and pressure due to the f(G) contribution defined as 

PG = Gfc - f - 24H 3 /g, (13) 

Pg = 16 H 3 f G + IGHHfc + 8H 2 fc - Gf G + /. (14) 


By using of Eqs. (13) and (14), one can obtain the equation of state (EoS) parameter due 
to the f(G) contribution as [14] 

Pg 16H 3 f G + 16HHf G + 8H 2 f G -Gf G + f 

WC = - = -:-• (15) 

Pg Gf G - f - 24ET 3 / g 

Also from Eqs. (9) and (10), the effective EoS parameter can be obtained as 


2 H 

UJeS- 1 3h2 - 

(16) 

Moreover, the continuity equations governing the pressureless matter (p m 
(p r = p r /3) and the f(G) contribution satisfy 

= 0), the radiation 

Pm T 3 H Pm = 0) 

(17) 

Pt + 4 Hp T = 0, 

(18) 

PG + 3 H (p G + pg) = 0. 

(19) 

Equations (7), (8), (17) and (18) determine the dynamics of the f(G )~gravity system (1) in a 


homogeneous and isotropic background. We will study this issue in section 4. 


3 GSL IN f(G )~GRAVITY 

Here, we are interested in exploring the GSL of gravitational thermodynamics in the context 
of f(G)- gravity. To this aim, we consider a spatially flat FRW universe filled with the matter 
and radiation. We further assume that the boundary of the universe to be enclosed by the 
dynamical apparent horizon with the Hawking temperature. For a spatially flat FRW universe, 
the dynamical apparent horizon takes the form [39, 40] 

r A = H- 1 , (20) 

which is same as the Hubble horizon. Following [24], the Hawking temperature on r A is given 
by 


T a 


l 

27rr A 


ta 

2 Hr a 


( 21 ) 


4 



Now we are going to use the first law of thermodynamics to find the general condition needed 
to hold the GSL in f{G)- gravity. The entropy of matter and radiation inside the horizon are 
given by the Gibbs equation [29] 


TA.d5 m — dF m + p m dV, 


( 22 ) 


TAdS' 1 - = dF r + p v dV, (23) 

where E m = p m V and E r = p v V. Also V = is the volume of the dynamical apparent 

horizon f a containing the pressureless matter (p m = 0) and radiation (p r = p r / 3). 

Taking time derivative of Eqs. (22) and (23) and using (17) and (18) one can get 


T a S = 4vn~A (Vm + ^Pr^j (f A - Hr A ) , 


(24) 


where S = 5 r + S m . Replacing p m and p r from Eqs. (7) and (8) into the above relation and 
using r\ = if -1 , one can obtain 


2vr 


T A S = 




H 4 


4 H - 16 H [h 2 - 2i/) / g + 1677 2 / g 


The horizon entropy in f(R, G)-gravity is given by [19] 

OR „ dG 


S A = -7T 


Fr 


OR, 




+ f g 


OR, 




^a/3^7 pdA, 


(25) 


(26) 


where Fr, = — , Fq = q^< G ' > and A = 47n'^ is the area of the apparent horizon (r A = 

H ~ 1 ). Also the quantity e a p is normalized as e a ^e a p = —2 and antisymmetric under the 
exchange a < — (5. It is the binormal vector to the bifurcation surface [19]. For the action (1) 
we have 

F(R,G) = R + 2f(G), (27) 

therefore Eq. (26) yields 

(28) 


5 a = 8tt 2 (H~ 2 + 8/ g ) • 


Taking time derivative of Eq. (28) and using the Hawking temperature (21), one can derive the 
evolution of the horizon entropy as 


T a S a = 


Summing up Eqs. (25) and (29), the GSL in f(G )-gravity yields 

T A 5tot = [2 H 2 + 8HH(dH + 3F 2 )/ g + 16F 2 (i7 + H 2 )f G 


(h + 2 H 2 ) 

/ • 0 • \ 


H 4 

(-2 H + 8H 3 f G ) . 

(29) 


(30) 


where 5 tot = 5 r + 5 m + S A - Equation (30) shows that in f(G) gravity, the validity of the GSL, 
i.e. FA5tot > 0, depends on the explicit form of the /(G) model. For the Einstein gravity 
(/(G) = 0), one can immediately find that the GSL (30) reduces to 


T A S tot = 


AttH 2 


>0, 


(31) 


which shows that the GSL is always fulfilled throughout history of the universe. In section 5, 
we examine the validity of the GSL (30) for two viable /(G)-gravity models. 
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4 Dynamics of f(G )~gravity 

To study the dynamics of a general f(G) model, we use the following dimensionless variables 
[14, 41] 


Gf G 
' ? i 3H 2 ’ 

(32) 

/ 

X ' 2 3i7 2 ’ 

(33) 

x 3 = -8 Hfe, 

(34) 

x . - Q - Pr 

4 r “ 3 H 2 ' 

(35) 

G H 

X5 ~~ 24i7 4 ~ H 2+ ' 

(36) 

x 6 = H. 

(37) 


With the help of above definitions and using Eqs. (7), (8), (17) and (18) one can get a set of 
first order differential equations governing a general f(G ) model as [14] 


dx’i 


— = -x 3 x 5 ( 1 •! — ) i 2.i:i (1 — .<- 5 ) 


m 


&X 2 X3X5 , , 

-TT7 = - + 2 X 2 (1 -X 5 ), 

div m 

= -^3(1 + x 5 ) + 2 x 5 + 1 - 3 (xi + x 2 ) + x 4 , 
d Jv 

^ = -2x 4 (l + x 5 ), 


dx 5 

dN 


= x 5 


4(1 - x 5 ) - 


X3X5 


mx 1J 


^ = (X 5 - 1 )X 6 , 

where N = ln(a/oj) and a t is the initial value of the scalar factor. Also 

G/gg 


m = 


/g 


(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

(44) 


Notice that the variable X6 = H, Eq. (37), and the differential equation (43) are absent in [14]. 
Using Eqs. (32)-(35), one can rewrite Eqs. (7) and (13) as 


U m = 1 — xi — x 2 — x 3 — X4, 

= xi + x 2 + x 3 , 

where fl m = and Qq = ■ Also we have Q m + O r + Qq = 1. 


(45) 

(46) 
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From Eqs. (15) and (16), one can rewrite cjg and 0J e fi in terms of the variables Xj as 


—2x5 — X4 — 1 

3(xi + x 2 + x 3 ) ’ 

1 , 

— ^(2X5 “1“ 1)* 

Also the deceleration parameter takes the form 


(47) 

(48) 


H 

q = ~ 1 ~Jp= 


(49) 


and 


H 6 f GG 


rnx 1 
192 X 5 ' 


(50) 


Notice that in the context of /(G)-gravity, the quantity H e f GG plays an important role. In a 
viable f(G ) model, the condition 0 < if 6 f GG < 1/384 is necessary in order to have a stable de 
Sitter point [15]. 

From Eqs. (13) and (14), the energy density and pressure due to the /(G) contribution can 
be expressed in terms of the variables x t as 


p G = 3 H 2 (xi + x 2 + x 3 ), 


(51) 


P G = -H 2 (2x 5 + x 4 + 1). (52) 

With the help of Eqs. (32)-(43), the GSL (30) reads 

T A S tot = 27r{2x 5 [3(xi + x 2 - 1) - (x 4 + x 5 )] + x 3 (5x 5 - 1) + 2 }. (53) 

Therefore for a given /(G) model, solving the set of first order differential equations (38)-(43) 
numerically, one can obtain the evolutionary behaviours of H, G, f2 m , Q G , cu G , w e g, q and 
r A Stot. In what follows, we investigate the dynamics of two viable /(G) models and examine 
the validity of the GSL, i.e. T A S tot > 0. 


5 Two viable f(G ) models 


Here, we are interested in investigating the GSL in two viable /(G) models. The first model has 
the form [14] 

/(G) = a(G 2 -$y , Model I, (54) 

where a and f3 are two constants of the model. The second /(G) model is given by [15] 


= arctan (^r'j ~ ^\/G* hi ^1 + ^j — aA\/Gh, Model II, 

where a is an arbitrary constant and A is a positive constant. Also G* = Hq and Hq is the 
Hubble parameter at present. 

With choice of suitable initial conditions, we numerically solve the differential equations 
(38)-(43) for both model I and model II. The evolutions of the Hubble parameter H , the Gauss- 
Bonnet curvature invariant term |G| and the quantity H 6 f G Q, Eqs. (36), (37) and (50), versus 
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N = ln(o/a;) for model I and model II are plotted in Figs. 1 and 2, respectively. Figures show 
that: (i) the Hubble parameter decreases during history of the universe, (ii) The Gauss-Bonnet 
curvature invariant term changes its sign when it transits from the standard radiation/matter 
dominated epochs to the accelerated era. (iii) The quantity H 6 /qq satisfies the condition 
0 < fGC,(Gi) < 1/384 which shows that both models have a stable de Sitter attractor, (iv) 
H , |G| and H 6 /qq at late times go to a constant value when the universe enters a de Sitter 
regime. Notice that the result of Fig. 2 for model II is the same as that obtained in [15]. 

The evolutions of the density parameters Q r , fi m , Hq and the effective EoS parameter w e g, 
Eqs. (45), (46) and (48), versus N for model I and model II are plotted in Figs. 3 and 4, 
respectively. Figures illustrate that: (i) for both models, H r , fl m , Hq and cj e ff behave like the 
ACDM model in the radiation/matter dominated epochs, (ii) For model I, w e ff oscillates rapidly 
during the accelerated epoch and goes deep into the phantom-like region as the universe enters 
the de Sitter period, (iii) For model II, cu e ff oscillates slowly around —1 as the system enters the 
epoch of cosmic acceleration, which implies that the de Sitter solution is a stable spiral. Note 
that the results of Figs. 3 and 4 are the same as those obtained in [14] and [15], respectively. 

The evolutionary behavior of the EoS parameter wq due to the f(G) contribution, Eq. (47), 
for model I and model II is plotted in Figs. 5 and 6, respectively. Figures present that cjg 
oscillates quickly at early times and approaches a de Sitter regime at late times, as expected. 

The evolution of the deceleration parameter q, Eq. (49), for model I and model II is plotted 
in Figs. 7 and 8, respectively. Figure 7 clears that for model I, the deceleration parameter starts 
from q = 1 corresponding to the radiation dominated epoch, then shows a cosmic deceleration 
(q > 0) to acceleration (q < 0) transition [42] and finally oscillates rapidly into the de Sitter 
regime (q = —1). Figure 8 presents that for model II, q varies from the matter dominated epoch 
(q = 0.5), then transits from a cosmic deceleration to acceleration and approaches smoothly a 
de Sitter regime at late times, as expected. 

In addition, we turn to check the energy conditions in both model I and model II. The energy 
conditions are as follows [43, 44, 45]: 

(i) The null energy condition (NEC), i.e. pq + pc > 0. 

(ii) The strong energy condition (SEC), i.e. pc + Pg > 0 and pc + 3pc > 0. 

(iii) The weak energy condition (WEC), i.e. pc + PG > 0 and pc > 0. 

(iv) The dominant energy condition (DEC), i.e. pc > 0 and pc > [pel- 

Using Eqs. (51) and (52) the evolutionary behaviors of pc + pg, PG + 3pg> PG and [pel 
versus N are plotted in Figs. 9—16. Figures 9, 11, 13 and 15 illustrate that for model I, the 
energy conditions during the standard radiation/matter dominated epochs are violated in some 
ranges of N. Thereafter, pc +Pg> Pg + 3pci PG and |pg| oscillate rapidly and finally the energy 
conditions, but the SEC, for model I are satisfied in the late times. Figures 10, 12, 14 and 16 
show that the energy conditions for model II behave like model I. But the difference is that 
in the future, Pg + Pg> PG + 3pc> PG and [pel for model II vary smoothly when the universe 
approaches a de Sitter regime. 

Finally, we examine the validity of the GSL for both models. In Figs. 17 and 18, we 

plot the variation of the GSL, Eq. (53), versus N for model I and model II, respectively. 

Figures illustrate that for both models, the GSL during the radiation/matter dominated epochs 
is fulfilled. Thereafter when the universe enters the cosmic acceleration era, i.e. q < 0 see Figs. 
7 and 8, the GSL does not hold (i.e. T\S to t < 0) in some ranges of N. At late times, the GSL 
for model I oscillates rapidly and for model II approaches smoothly into the de Sitter universe, 
adiabatically (i.e. Ta Stot = 0). 



6 Conclusions 


Here, we investigated the GSL of gravitational thermodynamics in the framework of /(G)- 

gravity. To do so, we considered a spatially flat FRW universe filled with the pressureless matter 

and radiation. We supposed the boundary of the universe to be enclosed by the dynamical 

apparent horizon with the Hawking radiation. We derived a general relation for the GSL which 

its validity depends on /(G)-model. Hence, for two viable /(G)-models containing /(G) = 
2 

a (g^ — p'j 3 [14] and /(G) = A-^g= arctan — ^VG* In (l + — aXy/G^ [15], we first 

solved numerically the set of differential equations governing the dynamics of /(G)-gravity. 
Consequently, we obtained the evolutionary behaviors of the Hubble parameter, the Gauss- 
Bonnet curvature invariant term, the pressureless matter, radiation and DE density parameters, 
the effective EoS parameter and the EoS parameter due to the /(G) contribution as well as the 
deceleration parameter. In addition, we turned to check the energy conditions containing the 
NEC, SEC, WEC and DEC. Finally, we examined the validity of the GSL for the two selected 
/(G)-models. Our results show the following. 

(i) The Hubble parameter H, the Gauss-Bonnet curvature invariant term |G| and the quan¬ 
tity H 6 f GG for both models at late times go to a constant value when the universe enters a de 
Sitter regime. Also both models have a stable de Sitter attractor, because H 6 f GG satisfies the 
condition 0 < H ( \ f GG (G\ ) < 1/384. 

(ii) The density parameters O r , Q m and Hg for both models behave quite similar to those of 
the ACDM model in the radiation and matter dominated epochs. 

(iii) The EoS parameters w e ff and uj g for both models start from the radiation/matter dom¬ 
inated epochs, then enter the phantom region (i.e. u < — 1) before reaching the de Sitter 
attractor with u = — 1 . 

(iv) The two selected /(G) models can give rise to a late time accelerated expansion phase 
of the universe. The deceleration parameter for both models shows a cosmic deceleration q > 0 
to acceleration q < 0 transition which is compatible with the observations [42]. Also for both 
models, q is ended with a stable de Sitter era (i.e. q —>• — 1). 

(v) The NEC, SEC, WEC and DEC for both models during the radiation/matter domi¬ 
nated eras are violated in some ranges of scale factor. But in the late times when the universe 
approaches a de Sitter regime, the all energy conditions, but the SEC, are satisfied. 

(vi) The GSL is fulfilled for both models during the standard radiation/matter dominated 
epochs. But when the universe becomes accelerating, the GSL is violated (i.e. Tj^Stot < 0) 
in some ranges of scale factor. At late times, the evolution of the GSL predicts an adiabatic 
behavior (i.e. TA<Stot = 0 ) f° r the accelerated expansion of the universe. 
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Model I 



Figure 1: The evolutions of the Hubble parameter H , the Gauss-Bonnet curvature invariant 
term |G| and the quantity H^/qq, Eqs. (36), (37) and (50), versus N = ln(a/ai) where a\ is 
the initial value of the scale factor. Auxiliary parameters are: a = and f 3 = —10^ 17 . 

Initial values are: x\ = —0.0025, X 2 = 0.005, X 3 = —0.01, X 4 = 0.99951, x§ = —0.99 [14] and 
x 6 = 7.95225 x 10" 7 . 



Figure 2: Same as Fig. 1 but for model II. Auxiliary parameters are: a = 10 and A = 0.075 
[15]. Initial values are: x\ = 189.249, X 2 = —189.248, X 3 = —0.0014, X 4 = 0.004, X 5 = —0.502 
and xq = 20 . 
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Model I 



Figure 3: The evolutions of fl m , flcj and w e ff, Eqs. (45), (46) and (48), versus N. Auxiliary 
parameters and initial values as in Fig. 1. 


Model II 



Figure 4: Same as Fig. 3 but for model II. Auxiliary parameters and initial values as in Fig. 2. 
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tion of the EoS parameter wq, Eq. (47), versus N for model I. Auxiliary 
rl values as in Fig. 1. 


Model II 



ig. 5 but for model II. Auxiliary parameters and initial values as in Fig. 2. 
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Model I 



Figure 9: The evolution of pq + po, Eqs. (51) and (52), versus N for model I. Auxiliary 
parameters and initial values as in Fig. 1. 


Model II 



Figure 10: Same as Fig. 9 but for model II. Auxiliary parameters and initial values as in Fig. 2. 
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Figure 15: The evolutions of pq and |pg|j Eqs. (51) and (52), versus N for model I. Auxiliary 
parameters and initial values as in Fig. 1. 


Model II 



Figure 16: Same as Fig. 15 but for model II. Auxiliary parameters and initial values as in Fig. 
2 . 


20 





































